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A solid is typically deemed amorphous, or glassy, when there are no Bragg peaks in its diffrac- 
tion pattern. We construct a two dimensional configuration of Ising spins with an autocorrelation 
function which vanishes at all nonzero distances, so that its scattering pattern is flat. This configu- 
ration is a ground state of a spin Hamiltonian with deterministic, translationally-invariant and local 
interactions. Despite ostensibly being amorphous, this configuration has perfect underlying order. 



The question of whether a system is considered ordered 
is an evolving subject. Before 1984, order was equated 
with crystallinity; thereafter, quasicrystals — aperiodic 
solids with perfect long-range translational order — were 
shown to be possible 1 . All other solids were regarded 
as glassy, or amorphous. It is natural to ask whether 
there are any other types of organization which may be 
perfectly ordered, while being neither crystalline nor qua- 
sicrystalline [3]. 

Recently, a new criterion for the definition and quan- 
tification of long-range spatial order was proposed [U [5]. 
This definition subsumes crystals and quasicrystals, and 
includes many other systems which would have been clas- 
sified as disordered due to their scattering spectrum. One 
example of this is the Rudin-Shapiro sequence [B], which 
may be thought of as a one-dimensional configuration 
of Ising spins, which has an absolutely continuous flat 
Fourier spectrum. Given that such configurations are 
mathematically possible, the next important question is 
whether they are physically relevant. That is: Can the 
configuration be the ground state of some reasonable 
Hamiltonian? 

In this Letter we develop a Wang tiling [7] which is 
a two-dimensional generalization of the Rudin-Shapiro 
sequence. One interpretation of this tiling is as a con- 
figuration of Ising spins (a = ±1) on a square lattice. 
Although its autocorrelation function is identically zero 
except at the origin, it is perfectly ordered in the sense of 
Reference [I]. Moreover, we shall show that this tiling is 
the ground state of a short-range Hamiltonian. Although 
this system may be thought of as somewhat artificial, the 
fact that it minimizes the energy resulting from a local 
interaction lends credence to the suggestion that physi- 
cal systems may realize ordered states which are neither 
crystalline nor quasicrystalline at low enough tempera- 
tures. Such states would not be distinguished by Bragg 
peaks in their scattering function, and thus might be in- 
correctly categorized as disordered. 

In order to develop the tiling [5], we will first consider 
several related polynomials [S]. Denoting the spin at the 
(square) lattice vertex fc, I by &k,i, we define the polyno- 
mial P(x, y) by 



where x,y are complex variables. Note that if x,y are 
defined on the unit circle, we may write x = e l<?x and 
y = e lqy , which, when substituted into P(x,y), yields 
the lattice Fourier transform of the {cfc,;}: 

^[M]=X>k,J e^e u ^ 
k,l 

Now, if the amplitude of the Fourier transform is a con- 
stant, independent of (q x ,Qy), the autocorrelation func- 
tion must vanish except at the origin: 

(Px,y&x+£,y+n) = S^fi^nfl (1) 

Thus the problem of finding a "deterministic amorphous" 
configuration of spins reduces to finding a deterministic 
polynomial in the variables x = e lqjl and y ~ e lqy which 
has a constant modulus. 

Let us now consider the polynomials P n = P n (x,y), 
and define them recursively using three auxiliary poly- 
nomials, Q n , R n , S n , where the dependence on x,y has 
been omitted for ease in reading: 

P n+ i = P n + x 2 Q n +y 2 R n + (xy) 2 S n 

Q n +i = P n + x 2 Q n -y 2 R n - (xy) 2 S n 

R n+ i = P n - x 2 Q n + y 2 R n - (xy) 2 S n (2) 

S n+1 = P n - x 2 Q n - y 2 R n + (xy) 2 S n 

Pq — Qo = Ro — So = i 

P n ,Q n , Rn, S n are all polynomials of degree 2 n+1 — 2 
and their coefficients take the values ±1. In particular, 
P n (x,y) = Yfk,i=o a ki ^V, where a kl = ±1. Thus we 
may consider a^i to be a configuration of Ising spins on 
a square lattice of linear size 2™, with the total number 
of spins being N — 4™. 

Assuming now that x and y have unit modulus, it is 
easily verified that 

+ R n+1 + S n+1 \ 2 = 16\P n \ 2 

\Pn+l + Qn+1 — Rn+1 — S n +l\ = l6\Q n \ 2 
\Pn+l — Qn+1 + Rn+1 — 'S'ri+l | 2 = 16|-R n | 2 
\Pn+l — Qn+1 — Rn+1 + S n +i\ 2 = 16\S n \ 2 

Summing these, we get 



k,l 



\Pn+l\ 2 + \Qn+l\ 2 + \Rn+l\ 2 + \S n +l\ 2 

= 4(|P„| 2 + |Q„| 2 + K| 2 + |S„| 2 ) 
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FIG. 1. Part of the configuration of Ising spins that is a gen- 
eralization of the Rudin-Shapiro sequence. Each black square 
corresponds to a spin a — 1, and each white square corre- 
sponds to a spin a — —1. Here < x, y < 64. 



FIG. 2. (color online) The configuration of spins, V, in 
terms of the eight lettters a+,a_, /?+,/?_, 7+, 7_,<5+,<5_ for 
< x, y < 256. Since each four spins are assigned one letter 
there are only 128 x 128 pixels in the figure. 



with |P | 2 + IQol 2 + l-Rol 2 H 
recursion whose solution is 



1 Sq 1 2 = 4. This is a simple 
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\Rr, 



\SJ 2 =4 r 



(3) 



This immediately implies that \P n \ < 2 n+1 , which scales 
as the square root of the number of spins; this implies 
an absolutely continuous Fourier spectrum [10], It is, in 
fact, possible to say more: in the thermodynamic limit 
the Fourier spectrum is flat, i.e. \P„\ is a constant inde- 
pendent of (fe, and 



IP, 



N, (as N -> 00) 
which are all of the form ±x ka,: 



(4) 




The terms in P n , which are all ot the torm ±x n "y", can 
be organized in the form of a 2" x 2™ matrix V n such 
that if u k jx k y l G P n (x,y) then (V n )l,k = a k,l- In this 
manner we also define matrices for the three auxilliary 
polynomials, and denote them Q n ,7?. n ,5 n . This allows 
us to write the recursive definition ^ as 



(5) 



with Vo — Qo = ^0 = S = 1. 

In each iteration, the linear size of the matrices 
Pn,Qn,'R-n,Sn is doubled. After n iterations each ma- 
trix has N = 4™ entries, which are ±1. We regard the 
matrix V n as a configuration of Ising spins on a square 
lattice; a configuration of 64 2 spins (n = 6) is shown in 

Fig.g 

Substitution Rules. Eq. ^ implies a recursive con- 
struction rule with which we may construct a Wang 
tiling. By looking at the matrix V n +i, we note that a 
configuration at generation n can be subdivided into four 



quadrants; each representing a configuration at genera- 
tion n — 1. These quadrants may again be subdivided. 
This iterative procedure "bottoms out" at generation, 
sav|12). n = 1. The configuration V n is now subdivided 
into 2x2 blocks of the form ±Vi,±Q 1 ,±TZ 1 ,±S 1 . If 
we simultaneously replace these 2x2 blocks by 4 x 4 
blocks as prescribed by Eq. ([5| we transform V n into 
V n +i- This procedure allows us to forgo configurations 
±Q n ,±lZ n ,±S n for n > 1 and inflate the configuration 
V n -> V n +i directly. 

Let us now denote each of the eight 2x2 blocks by a 
letter from the alphabet {a + , a_, /3_, 7+, 7_, <5 + , <5_} 
as follows 



a + = 
P+ = 

7+ = 



-a- = Vi 
-P- = Ci 

-7- - ^1 

: -5- = 5l 
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(6) 



Each of these letters will represent a set of tiles in the 
Wang tiling. In terms of these letters, the configuration 
V n is shown for n — 8 in Figure [2] For ease in visu- 
alization, we will represent each of the tiles as a 2 x 2 
'semaphore' block, with 1 represented by a black square 
and -1 by a white square, as shown in Figure [3] 

By Eq. ([5]), each of the above 2x2 blocks is to be 
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FIG. 3. The eight basic tiles in the 'semaphore' representa- 
tion, in which a black square represents a spin a = 1, derived 
from Eq. B. 



if the markings on their shared edge match — these are the 
so-called matching rules. We note that although the con- 
figuration above is generated substitutionally from eight 
basic tiles, there will be more than eight Wang tiles in 
terms of their markings. 

By associating an energy penalty +J for each edge 
with mismatched markings, we define a Hamiltonian 



H = j£> 



S(T„ 
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substituted by a 4 x 4 block in the following manner 



(7) 

ct-t- 0- a_ Ox. 

7+ 



These substitution rules are shown in Figure [4] The sub- 
stitution rules obey two important properties: (i) each 
of the pairs (aq_,a_), (/?+,/?_), (7+, 7-), and (<5 + ,<5_) 
occupy a specific quadrant on the rhs of each substitu- 
tion rule, (ii) in each pair, one partner is the negative of 
the other, and its substitution rule is the negative of its 
partner's substitution rule. By repeatedly applying these 
rules, we can generate a tiling of arbitrary size. The am- 
plitude of the Fourier transform of configurations arising 
from substitution rules obeying these two properties was 
rigorously shown [TT] to be flat, independent of wave vec- 
tor. 
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FIG. 4. The two-dimensional substitution rules, correspond- 
ing to Eq. |7f . These substitutions generate the configuration 
V when iterated on a 2 x 2 block of all black squares. 



Tiling model. By a theorem due to Mozes [13] , config- 
urations emerging from substitution rules such as those 
111 Eq. may be enforced by a set of Wang tiles — unit 
squares with marked edges. The plane is to be tiled by 
copies of these tiles 14J such that two tiles may abut only 



where 5(0) = 1 and where T is a tile with edge markings 
T N , T s , T E , T w . The ground state of this Hamiltonian 
has zero energy and corresponds to a tiling in which there 
are no mismatches. 

We note that this system may be regarded as a gen- 
eralized Potts model. Let t be a Potts spin, taking M 
possible values, where M is the number of distinct Wang 
tiles (by Mozes' theorem, M is finite). We denote these 
values by t = 1 . . . M. Since two nearest-neighbor tiles 
have an interaction energy (+ J or 0), the Hamiltonian ffify 
can be written as 



+1,2/ ; 



(9) 



where Q N and Q EW are two M x M matrices whose 
entries are and 1. In a regular Potts model each Q 
matrix would have zeros on the diagonal and ones off 
the diagonal; here the matrices reflect the matching rules 
between pairs of tiles in their various juxtapositions. 

In order to fully define the system, we must construct 
the matrices Q used in Eq. |9]) , or equivalently, construct 
a set of Wang tiles which allow our configuration as the 
only legal tiling of the plane. We have constructed such 
a set by the two-step procedure described below. 

Step 1. Apply Mozes' construction [T3j to our 
tiling. In his construction there are two funda- 
mental types of tiles. The first type corresponds 
to the substitutional ("letter") tiles (in our case 



a + ,tv_,/3 + ,^_,7 + ,7_,(5 + ,(5_). 



These tiles are con- 



strained to occupy the sites of a diluted lattice of double 
the spacing of the physical lattice (e.g., sites where x and 
y are both even). The second type corresponds to tiles 
whose function is to synchronize [13] the letter tiles so 
that those will order into the tiling shown in Fig. [2] The 
"synchronizer" tiles are not allowed to occupy the sites 
of the diluted lattice. This construction leads a tile set 
consisting of 116 tiles [TS], of which 8 are letter tiles and 
108 are synchronizer tiles. The Hamiltonian is not yet 
that of Eq. ^ because different tiles are restricted to 
different sites of the lattice. 

Step 2. Define "supertiles" , each a plaquette con- 
sisting of one letter tile and the three synchronizer tiles 
lying to its right, bottom, and bottom-right. The edge- 
markings of the supertiles are the superposition of the 
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edge-markings of its constituent tiles. A careful com- 
binatorical counting shows that this step generates 400 
supertiles, for which the Hamiltonian is that of Eq. (J8|. 

The 400 tiles differ by their edge mark- 
ings but represent variations of the letter tiles 
{oq_,a_,/3 + ,/3_,7 + ,7_,5 + ,<S_}. By ^j, the tiles 
{a+, /3+, 7+, 5 + } are the negatives of the tiles 
{«-, /9_, 7-, 6—}, and indeed if we color the first 
set black and the second set white, the ground state of 
the Hamiltonian (|8| is the configuration shown in Fig. [TJ 

Underlying Order. Although its two-point function 
vanishes (JTJ), our tiling is deterministic, and has subex- 
tensive patch entropy [4]. Thus, we expect there to be 
some measure which would differentiate between it and a 
random arrangement of the same tiles. It is important to 
note that there is no a priori way of knowing what mea- 
sure would reveal the order of an arbitrary given tiling, 
but, for the system we examine in this paper, we have 
identified one such four-point correlation function 

n(k, I) — (7k,l<Tk+i,l<Tk,l+iO'k+i,l+i 

which is useful in revealing the underlying hierarchical 
order of our tiling. 

This is essentially a parity function as it returns +1 
when there are an even number of a = —1 spins in 
the quadruplet of spins in which k,l is the upper left 
position, and —1 when there are an odd number of 
a = — 1 spins. This function is plotted in Fig. [5] The 
average of this four-point correlation function is (tt) = 
(1/N) ^2 k l ir(k, I), and by direct inspection of Fig. [5l we 
see that one third of the sites of the the ground state nave 
7T = 1 and the other two thirds have ix = — 1. Therefore 

(tt) = - 1 
V 1 3 

and the four-point correlations, unlike the two-point cor- 
relation, are not generically zero. As may be seen from 
Fig. [5j the function ir(k, I) may be regarded as a super- 
position of lattices. The first lattice has a two-site spac- 
ing between sites with tt — 1. The second lattice has 
a four-site spacing between sites with n = 1, etc. This 
is a superposition of lattices in which the level n lattice 
(n > 1) has a lattice constant of 2™. 

We have constructed a deterministic two-dimensional 
system of Ising spins whose ground state is glassy in the 
sense that its autocorrelation function vanishes exactly 
in the thermodynamic limit. This configuration may be 
generated by a recursive method, which guarantees that 
its patch entropy [J] scales logarithmically with patch 
size. In this sense, the ground state configuration is or- 
dered. A tiling which enforces a hierarchical structure, 
albeit on a triangular lattice, has recently been studied 
by Byington and Socolar |16j and shown to exhibit an 
infinite sequence of phase transitions; The finite temper- 
ature properties of our system should also prove intrigu- 
ing, and will be detailed in a future publication. 




FIG. 5. The parity function ir(k,l) for < k, I < 16. Black 
squares represent 7r(k,l) = 1 and light gray squares represent 
7r(fc,0 = -1. 
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